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Abstract 

The correct form of the Schmidt decomposition of the stationary wave functions for a 
system of two interacting particles trapped in a two-dimensional harmonic potential is given. 



Keywords: Schmidt decomposition 



We note that in Ref.(H the expansion of the two-particle wave function (Eq.ll) is mis- 
takenly interpreted as the Schmidt decomposition. The mode functions appearing in it 
cannot be considered as the Schmidt orbitals due to their incorrect normalization in radial 
direction in two-dimensional (2D) space. It should be stressed that this mistake does not 
affect the validity of the results presented in Ref. [l[ since Eq.ll has not been used to their 
determination. However, because of the utility of the model in many areas of physics it is 
important to provide the correct Schmidt form decomposition that we derive below. 

The system of two interacting particles trapped in a 2D isotropic harmonic potential, 
irrespectively of the interaction potential between the particles, possess the stationary wave- 
functions (m r = M c = 0)0 that depend only on p and g, where \g\ = \{pi + p%)/2\ = 
q(Pi, P2,cos(ip 2 - (fi)) and \p\ — |p 2 — pi\ — p(pi, p 2 , cos(ip 2 - <pi)). Being the function of p 
and g only, ^ is symmetric under permutation of the particles ((pi,<pi) < — > (^2,^2)) and 
may be assumed to be real for simplicity. An application the method of Ref. Q to this case 
results in 

*(p, e ) = ]T Am(Pl,P2) e im ^e-^ 2 , (1) 

m=— 00...00 

where to ensure correct normalization in radial directions we have introduced ^p\p2 which 
is the crucial difference from Eq. 11 of Ref. [lj. The function A m (pi,p 2 ) is given by the 
following integral 

2tt r 2n 



lr The Hamiltonian of two interacting particles in the isotropic harmonic trap is separated into cen- 
ter of mass (cm.) and relative (rel) motion. The stationary states of this system may be chosen as 
^n,m r ,N,Mc(PuP2) = ^ eL (p)e im ^'- ip% m - (g)e iM ^ m - , where the functions and ip% m - are solutions of 
the radial Schrodinger equations rel. and cm., respectively. 
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where the simplification has been achieved by elementary symmetry considerations. Being 
real and symmetric, the function A m (pi, p 2 ) has the following Schmidt form 

A m (Pi,P2) = E^-xi m) (Pi)xi m) (p 2 ), (3) 
where the coefficients K s>m and the orbitals \s (p) satisfy the integral equation 

A m {pi,P2)x ( s n) {p2)d P 2 = ^, m xi m) (Pi)- 



The family {x s m \p)}s=o forms a complete and orthogonal set (f£° xi^X^dp = C5 SS >). 
Using the expansion ([3]) the final form of the decomposition of the wave function \1/ now 
reads 

(pl,¥l)v*(P2,(f2), (4) 
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where 

VsApM^^-t 1 ^- (5) 
VP 

Since the orbitals ([5]) satisfy the condition of orthogonality in 2D space 

roc r2-rr roc , . 

/ / PKrn V s',m' d P d ¥ = 2n5 ' Xs X? d P = 27rC5 mm' 5 ss' , 
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we recognize them as the Schmidt modes (natural orbitals). One can point out that they 
are the eigenf unctions of the angular momentum operator L z . 

In conclusion, we have presented in details a procedure of obtaining the Schmidt decom- 
position of the two-particle wave function that is a function of distances p and g only. 
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